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ABSTRACT 



t <~| Aims. In this study we try to answer the question whether or not subluminous Type la supernovae have additional distinctive properties when 
Oh examined from the point of view of the explosion geometry. 

^ Methods. We have performed optical spectropolarimetric observations of the Type la SN 2005ke at 3 epochs (days -8, -7, and +76). The 
< I explosion properties are derived by comparing the data to explosion and radiation transfer models. 

Results. The SN shows polarimetric properties that are very similar to the only other subluminous event for which spectropolarimetry is 
C^j available, i.e. SN 1999by. The data present a very marked dominant axis, which is shared by both the continuum and lines such as Sin ,46355, 
1 1 suggesting that the relatively large, global asymmetry is common to the photosphere and the line-forming region. The maximum polarization 
degree observed in the Sin /i6355 absorption reaches 0.39±0.08%. At variance with what is seen in core-normal Type la, SN 2005ke displays 
^ significant continuum polarization, which grows from the blue to the red and peaks at about 7000 A, reaching ~0.7%. The properties of the 
qq polarization and flux spectra can be understood within the framework of a subluminous delayed-detonation (DD), or pulsating DD scenario, 
or WD mergers. The difference in appearance with respect to core-normal SNe la is caused by low photospheric temperatures in combination 
ryQ w i m layers of unburned C, and more massive layers of the products of explosive C and O burning. The comparatively high level of continuum 
i— H polarization is explained in terms of a significant global asymmetry (~15%), which is well reproduced by an oblate ellipsoidal geometry within 

the general context of a delayed-detonation explosion. 
f—.j Conclusions. Our results suggest that SN 2005ke arose either from a single-degenerate system in which the WD is especially rapidly rotating, 
close to the break-up velocity, or from a double-degenerate merger. Based on the current polarization data, we cannot distinguish between these 
t—{ two possibilities. Possible tests are discussed. 



. £^ Key words, supernovae: general - supernovae: individual: 2005ke - ISM: dust, extinction 

X 

H 

1. Introduction These objects, commonly dubbed "91bg-like" after the pro- 
totypical event SN 1991bg (Filippenko et al. |1992| Leibundgut 
Subluminous Type la Supernovae (SNe la) are still a puzzle in et al. 1 1993b . deviate quite significantly from the behav- 
our understanding of this important class of explosive events ior denned by core-normal events (Branch et al. 2006). 
(Hillebrandt & Niemeyer|2000|l. Althoug h they contribute to Photometrically, they are about 2 magnitudes fainter, display 
"only" 15% of all Type la events (Li et al. 12011) , understand- rapidly declining light curves with narrower peaks, lack a see- 
ing their progenitors and the mechanism that powers their ex- ondary maximum in the NIR bands, have redder colors at 
plosions is fundamental to fully grasp the T ype la supernova maximum light, and obey a different relation between maxi- 
phenomenon (see, for instance, Pakmor et al. 12010) . mum luminosity and light curve shape (Garnavich et al. 12004) . 

Spectroscopically, they show cool spectra, characterized by 

intermediate-mass elements, particularly strong Oi and Tin 

Send offprint requests to: F. Patat absorptions, and lower expansion velocities. In addition, both 

* Based on observations made with ESO Telescopes at the Paranal iron-group elements and silicon are spread over velocities span- 
Observatory under Program IDs 076.D-0177 and 076.D-0178. 
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Fig. 1. A DSS image of NGC 1371. The location of SN 2005ke 
is marked by the cross. The arrow indicates the ISP position 
angle (0/s/>=17O degrees; see Sec.|4|. 



ning a large fraction of the ejecta, implying very substantial 
mixing (Taubenberger et al. 2008), as opposed to what hap- 
pens in core-normal objects. The estimated mass of 56 Ni is as 
low as 0.1 M Q (Stritzinger et al. 120061 1 . 

All these facts pose a challenge and call for a different 
explosion mechanism for subluminous SNe la than for core- 
normal SNe la. As of today, although several theoretical sce- 
narios have been put forward (see Pakmor et al. 120111 for a 
recent review), we do not have a clear idea about the origin 
of this sub-class of Type la SNe. On the observational side, 
there are indications that these objects form a rather homoge- 
neous and distinct class; however, there is growing evidence 
that transition objects, sharing properties of both subluminous 
and core-normal SNe, d exist (Maguire et al. 1201 11 1. 

An additional channel for investigating the diversity of 
SNe is polarimetry, which enables the study of the explosion 
geometry (Wang & Wheeler 120081 . Well studied 91bg-like 
events are quite rare, both because of their lower luminosity 
and their intrinsic lower rate. For the same reasons, sublumi- 
nous events studied through spectropolarimetry are exception- 
ally rare (Wang & Wheeler 2008). More precisely, so far the 
geometrical properties of the explosion have only been studied 
for one such SN, i.e. SN1999by, which showed very distinct 
features (Howell et al. |200TT ). 

In this paper we present a spectropolarimetric study of an- 
other 91bg-like event, SN 2005ke. This SN was discovered by 
Baek, Prasad & Li (120051 ) on 13 November 2005 (UT) in the 
Virgo Cluster spiral galaxy NGC 1371 (see Fig.[T|, located at a 
distance of 25.8+3.2 Mpc (Wood-Vasey et al. 120081) and reced- 
ing with a velocity Vi, fl /=1463 km s (Koribalski et al.|2004). 
A few days later, the object was classified as an underluminous 
Type la event caught before maximum light (Patat et al. 2005 ). 
Spectroscopic and photometric data were presented by Immler 
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Fig. 2. Spectroscopic evolution of SN 2005ke during the 
epochs covered by the observations presented in this paper. 
Phases refer to B maximum light (JD=2,453,699.16; Folatelli 
et al. 2010). For presentation the spectrum of the second epoch 
is shifted by logF=-0.5. 

Table 1. Log of VLT-FORS1 spectropolarimetric observations 
of SN 2005ke. 



Date 


MJD phase 


V 


airmass 


exp. time 


(UT) 


(JD-2,400,000.5) (days) 


(*) 


(average) 


(seconds) 


2005-11-16 


53690.10 -8.6 


15.2 


1.1 


4x300 


2005-11-17 


53691.10 -7.6 


15.0 


1.2 


4x (300+600) 


2006-02-08 


53775.06 +76 


17.7 


1.2 


4x900 



(*) V magnitudes are derived from Contreras et al. d20 lOt 



et al. d2006l Bufano et al. d2009l , Folatelli et al. d20T0b . and 
Contreras et al. (120101 ). to which we refer the reader for the 
general properties of this object. 

The paper is organized as follows. In Sect.|2]we present and 
discuss our observations and data reduction. The spectroscopic 
evolution during the phases covered by our observations is an- 
alyzed in Sect. [3] while Sect.[4]presents the spectropolarimetric 
data. The data are compared to explosion and radiation transfer 
models in Sect. [5] In Sect. [6] we discuss our results, while in 
Sect. [7] we summarize the conclusions. 

2. Observations and Data Reduction 

We have observed SN 2005ke on 3 different epochs, us- 
ing the FOcal Reducer/low-dispersion Spectrograph (hereafter 
FORS1), mounted at the Cassegrain focus of the ESO-Kueyen 
8.2m telescope (Appenzeller et al. 1998 ). In this multi-mode in- 
strument, equipped with a 2048x2048 pixel (px) TK2048EB4- 
1 backside-thinned CCD, polarimetry is performed introducing 
a Wollaston prism (19" throw) and a super-achromatic half- 
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Fig. 3. Spectroscopic comparison of subluminous SNe 2005ke, 
2005bl (Taubenberger et al. 120081 and core-normal SN 1994D 
(Patat et al. 120061) at -7/-6 days. 



wave plate. In order to reduce some known instrumental prob- 
lems (see Patat & Romaniello 2006) we have always used 4 
half-wave plate (HWP) angles (0, 22.5, 45 and 67.5 degrees). 
The SN was observed at two pre-maximum phases (-8.6 and 
-7.6 days) and one additional epoch (two and a half months 
past maximum light). Exposure times ranged from 5 to 15 min- 
utes per plate angle and, on the second epoch, the sequence 
was repeated in order to increase the signal-to-noise ratio. All 
spectra were obtained with the low-resolution G300V grism 
coupled to a 1.1 arcsec slit, giving a spectral range 3300-8600 
A, a dispersion of ~2.9 A pixel -1 and a resolution of 12.4 
A (FWHM) at 5800 A. Data were bias, flat-field corrected 
and wavelength calibrated by means of standard tasks within 
The RMS error on the wavelength calibration is about 
0.7 A. The wavelength scale was corrected to the rest-frame us- 
ing the host galaxy recession velocity (1463 km s _1 ; Koribalski 
et al. 120041 . 

The ordinary and extraordinary beams were processed sep- 
arately. Stokes parameters were computed by means of spe- 
cific routines written by us, and error estimates were performed 
following the prescriptions described by Patat & Romaniello 
(2006), while the HWP zeropoint angle chromatism was cor- 
rected using tabulated data (Jehin, O'Brien & Szeifert 2002). 
In order to increase the signal-to-noise ratio, multiple data sets 
obtained at the same epoch were combined, and the final Stokes 
parameters binned in ~52 A wide bins (20 pixels). 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Science 
Foundation. 
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Fig. 4. Evolution of Sin /16355 for a compilation of underlu- 
minous Type la: 1991bg (Turatto et al. 119951 . 1998de (Modjaz 
et al. l200Tl Matheson et al.|2008 ), 1999by (Vinko et al. l200"Tl 
Howell et al. [2Q0Tb . and 2005bl (Taubenberger et al. 120081 . The 
dashed line is a best fit to the data after excluding SN1998de. 

Flux calibration was achieved through the observation of 
spectrophotometric standard stars with the full polarimetric op- 
tics inserted (HWP angle set to degrees). Instrumental polar- 
ization and the position angle offset were checked by observing 
polarized and unpolarized standard stars, obtained within the 
FORS 1 calibration plan. 

The log of observations is reported in Table [T| where the 
SN phases were computed with respect to the B maximum light 
(JD=2,453,699.16; Folatelli et al. l20T0l ). 

3. General properties of SN 2005ke 

From the classification spectrum it was clear that SN 2005ke 
was a subluminous Type la event (Patat et al. 2005 ), very sim- 
ilar to SN 1999by (Howell et al. llXJTl Garnavich et al. 120041 . 
The characteristic Ti n feature at about 4200 A was clearly de- 
tected in the pre-maximum spectra. The spectroscopic evolu- 
tion of SN 2005ke during the epochs discussed in this paper 
is presented in Fig. [2] The close resemblance to another well- 
studied subluminous SN la 2005bl (Taubenberger et al. 2008 ) 
is illustrated in Fig. [3] where for comparison we have also plot- 
ted the spectrum of the core-normal SN 1994D at the same 
pre-maximum epoch. The expansion velocity, deduced from 
the minimum of the Si n /16355 absorption trough on day -7, 
is 11,300 km s _1 , very similar to the values measured for the 
subluminous SN 2005bl (Taubenberger et al. 120081 . On the 
following day, the velocity dropped to about 10,900 km s _1 . 
These values are compared to those published for other sublu- 
minous la events in Fig. |4] Although only early measurements 
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are available for SN 2005ke, and therefore no firm conclusion 
can be drawn, this object seems to conform to the behavior 
shown by the other 91 bg-like events. A best fit to the combined 
data of SN 1991bg, 1999by and 2005ke gives a velocity gradi- 
ent v=106±9kms~ 1 day -1 . 

SN 1998de displays higher velocities in the pre-maximum 
epochs, while a month past maximum light it falls very close to 
the relation derived from the previous objects. Its average ve- 
locity gradient is ~ 175 km s _1 day -1 . Although this indicates a 
possible break in the homogeneity otherwise seen in this sub- 
class of SNe la (Taubenberger et al.[2008), all the objects plot- 
ted in Fig. ^definitely qualify for the HVG class introduced by 
Benetti et al. d20051 >. 

From a photometric point of view, SN 2005ke has shown 
the same properties as other 91bg-like objects (see Contreras 
et al. 120101 1, i.e. the absence of a secondary peak in the 
red passbands, an intrinsically red color and a fast decline 
rate Awi 5 (fi)=1.76±0.01 (Folatelli et al. EOTOl . For compari- 
son, SN 1999by and SN 2005bl reached Am 15 (B)= 1.90+0.05 
(Garnavich et al. |2004] | and 1.93+0.10 (Taubenberger et al. 
2008), respectively. SN 2005ke reached a peak magnitude of 
V-14.2 (Contreras et al. 2010), which coupled to the esti- 
mated distance modulus yU=32.06+0.27 (Wood-Vasey et al. 
120081 yields an absolute magnitude My=-17.9+0.3. We note 
that a fainter magnitude is obtained (My ~ -17 ) if one uses 
the Tully-Fisher distance modulus (ju=3 1.16, Tully[l988). 

4. Spectropolarimetry 

Although we obtained two separate pre-maximum epochs 
(days -8.6 and -7.6), a close analysis of the data shows that 
there is no statistically significant evolution between the two 
dates. Therefore, in view of the subluminous nature of the ob- 
ject, we have merged the two data sets in order to maximize 
the signal-to-noise ratio. The result is shown in Fig. [5] where 
we plot the combined spectropolarimetric data set (correspond- 
ing to a total exposure time of 1200 s per half-wave plate angle) 
on the Q - U plane. 

The data show a dominant axis (Wang et al. 1200 It Wang 
& Wheeler [20081 l. not passing through the origin of the Q - U 
plane, hence indicating the presence of different components, 
characterized by different polarization degree and position an- 
gle (see for instance Patat et al. [2009b, their Appendix B). The 
polarization variation along the dominant direction appears to 
be moderate, but significant, reaching about 0.8%. In contrast, 
the polarization variation in the perpendicular direction is con- 
sistent with the measurement errors (see Sect. |4.2] >. 

The first step for disentangling the different contributions is 
the subtraction of any interstellar polarization (ISP) that might 
arise within aligned asymmetric dust grains along the line of 
sight. 

4.1. Reddening and interstellar polarization 

The Galactic reddening along the line of sight to SN 2005ke 
is £Vi/=0.023 (Schlegel, Finkbeiner & Davis [19981 . Using 
the relation found by Serkowski, Matheson & Ford ( 1975) this 



Table 2. ISP determinations 



Range 


phase 


ID 


(Qu P ) 


(U isp ) 


(A) 






(%) 


(%) 


4800-5600 


pre-max 


A 


+0.007±0.012 


-0.186±0.017 


4000-4200 


pre-max 


B 


+0.155±0.040 


-0.053+0.047 


5100-5300 


pre-max 


C 


-0.029±0.025 


-0.303+0.030 


4500-7500 


late 


D 


+0.143±0.029 


-0.054+0.029 



implies that the Galactic interstellar polarization (ISP) is ex- 
pected to be lower than ~0.2%. The compilation of polarimet- 
ric data for Galactic stars by Heiles (1999) contains only one 
entry within 3 degrees from NGC 13710 i.e. HD 22332. This 
star, at a projected distance of 1.1 degrees from SN 2005ke, 
has a linear polarization P=0.13% at position angle PA=4.1 
degrees (Matheson & Ford 1970), in good agreement with the 
low Galactic extinction along this line of sight. Since the ex- 
tinction within the host galaxy is estimated to be comparably 
low (£ B _ V =0.036±0.005; Folatelli et al. EOTOl l. the total ISP is 
expected to be of order of a few 0.1%. 

Several methods to derive the level of ISP in SN spectra 
have been proposed in the literature (see for instance Leonard 
et al. [20001 Howell et al. |200T1 Wang et al. |200H and used 
in a number of cases in the SN field. They are based either 
on the analysis of the Q - U plane, or on the assumption that 
some portions of the SN spectra are intrinsically unpolarized 
(Wang, Wheeler & Hoflich [T997l Howell et al. l20%n Chornock 
et al. 2006), so that a non-null continuum polarization can be 
interpreted as arising within the interstellar medium. 

Although the application of these methods is possible for 
this object (see below), we attempted an alternative solution, 
which is trying to actually measure the ISP. For this purpose we 
have obtained spectropolarimetry of SN 2005ke at a rather late 
epoch, about two and a half months past maximum light, when 
the SN is entering the nebular phase. Since at these epochs the 
continuum polarization from an asymmetric photosphere is ex- 
pected to be null, any residual continuum polarization can be 
interpreted as arising in intervening interstellar dust. The very 
well-studied case of SN 2006X (Patat et al. 2009a) showed that, 
with only the remarkable exception of the Ca n triplet, line po- 
larization is null at epochs later than one month past maximum 
light. This allows one to use a wide wavelength range in order 
to reduce the uncertainty. 

The spectropolarimetry of SN 2005ke on day +76 is pre- 
sented in Fig. [6] The data do not show any statistically sig- 
nificant polarization variation across the whole spectral range. 
The weighted average Stokes parameters, estimated in the 
wavelength range 4500-7500 A, are (Qi S p}=0. 14+0.03% and 
(I//sp)=— 0.05±0.03%, corresponding to a polarization degree 
Pisp=®- 15+0.03%, and a position angle 8 IS p=110±6 degrees. 
The rms deviation of single data points from the weighted av- 
erage value is 0.26%, i.e. fully compatible with the measure- 
ment errors. Because of the weakness of the SN at this epoch 
(V ~ 17.7, i.e. almost 3 magnitudes fainter than at maximum), it 
is not possible to deduce the wavelength dependency of the ISP. 
For a Serkowski- Whittet law (Serkowski et al. 119751 Whittet 

2 http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=II/226 
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Fig. 5. Combined spectropolarimetry of SN 2005ke (November 16 and 17, 2005). The large filled circle marks the ISP estimated 
from the late-time data (see Sect. 



4.1 



while the circled crosses indicate B, A and C ISP estimates (from right to left respectively; 
see Table [2}. The dashed line traces the best fit dominant axis, while the ellipse indicates the axial ratio and orientation of the 
data (see Sect. 4.2 1. For presentation, the size of the ellipse was arbitrarily scaled. The dotted circles indicate polarization levels 
of 0.25, 0.50 and 0.75%. 



et al. 1 1992| l the polarization variation is very mild across the 
optical domain, and its absolute amount is proportional to the 
maximum polarization. For the present case, this implies a very 
small (<0.05%) peak-to-peak variation across the spectral do- 
main covered by our observations, certainly below the accu- 
racy one can achieve with instruments like FORS (Patat & 
Romaniello 120061 . 

In addition, we note that the ISP wavelength dependency 
observed in several supernovae clearly deviates from the 
Galactic Serkowski law (Leonard & Filippenko 2001 1 Leonard 
et al. 120021 Maund et al. I2007al Patat et al. |2009al |2009bl 
Maund et al. 2010a). Therefore, rather than using a Serkowski 
law to model the ISP, we ignored its wavelength dependency 
and vectorially subtracted the constant values derived above. 
Once placed on the Q - U plane, the derived ISP happens to 
be very close to the dominant axis derived fitting all the data of 
the combined spectropolarimetry obtained in the pre-maximum 
epoch (Fig. [5] large filled circle). We also notice that the mea- 
sured ISP is very close to one of the two possible solutions one 
would consider in the absence of late-time data, at either end 
of the data along the dominant axis (see, for instance, the ISP 
derivation for SN 1999by, Howell et al. 1200 11 1. More precisely, 



the measured ISP essentially coincides with the solution imply- 
ing the lowest ISP, which is fully in line with the low reddening 
suffered by the SN. 

As an independent check, we have derived the ISP measur- 
ing the polarization in the wavelength ranges 4800-5600 A (A), 
4000-4200 A (B), and 5100-5300 A (C) on the combined pre- 
maximum spectrum. These ranges were proposed by Howell 
et al. (GOOD, Chornock et al. d2006b and Maund et al. d2010al > 
respectively, based on the assumption that in these regions the 
line blanketing is sufficient to produce complete depolarization. 
The results are compared to those obtained from the late-time 
spectrum in Tableland Fig. [5] (circled crosses). The ISP de- 
rived from the three wavelength ranges is <0.3%, which is con- 
sistent with the low level of polarization expected from the mild 
extinction along the line of sight to SN 2005ke. The estimate 
obtained for the B range coincides (within the statistical errors) 
with the late-time measurement, indicating that in this particu- 
lar spectral region the depolarization assumption holds. 

This is not the case (at least for SN 2005ke) for the other 
two wavelength intervals, which give ISP values that are in- 
consistent with the late-time data. In addition, the Stokes pa- 
rameters derived from range A place the ISP very close to the 
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5000 6000 7000 

Restframe Wavelength (A) 



8000 



Fig. 7. Decomposition of polarization in the combined pre-maximum data set along the dominant (central panel) and orthogonal 
(lower panel) axis. The upper panel shows the un-binned flux spectrum and main line identifications; the thin curve traces 
a synthetic telluric absorption spectrum. The horizontal dotted lines in the lower panel indicate the rms deviation from null 
polarization. 



centroid of the data set (see next section), so that notwithstand- 
ing the clear existence of a dominant axis, the implied polar- 
ization angle would become erratic (this would also turn into a 
very low intrinsic polarization level). As for range C, it must be 
noticed that this was proposed for the peculiar SN la 2005hk, 
where this particular region was indeed expected to have null 
intrinsic polarization (Maund et al. [201 0a). This clearly does 
not apply to SN 2005ke. 

We will adopt the ISP derived from the late-time data 
throughout the paper. 

Since dust grains tend to be aligned along the direction of 
the galactic magnetic field, which in turn follows the spiral pat- 
tern (see Scarrot, Ward-Thompson & Warren-Smith 1987), the 
ISP is expected to be tangential to the local spiral arm of the 
host galaxy. This has been observed, for instance, in SN 2001el 
(Wang et al. 120031 SN 2006X (Patat et al. I2009al), i n the 
bright optical transient in NGC 300 (Patat et al. |2009bl l, and 
in SN 2005hk (Maund et al. l2010ab . When the ISP within the 
host galaxy and the Milky Way are comparable, as is the case 



Table 3. Best fit dominant axis parameters 



Least Squares Fit Analysis 



a 


-0.19% (0.01%) p 


1.10(0.13) 


X 2 /DOF 


140/85 r„ 


0.43 


a 


0.16% e d 


113°.9(1°.6) 


Principal Components Analysis 


(Q) 


0.00 (0.01) <{/) 


-0.20 (0.01) 


b/a 


0.47 e d 


114.0 



Note: values in parentheses indicate rms uncertainties. 



for SN 2005ke, the position angle of the resulting polarization 
vector can differ substantially from the direction defined by the 
local spiral pattern. This is the likely explaination why the de- 
rived ISP position angle is far from being parallel to the spiral 
arm of NGC 1371 at the location of SN 2005ke (Fig. [T}. The 
ISP position angle (170+6 degrees) is not too different from 
that reported for the nearby star HD 22332 (4.1 degrees; see 
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4000 5000 6000 7000 
Wavelength (A) 



8000 



Fig. 6. Stokes parameters Q (central panel) and U (lower 
panel) on day +76. For presentation, the data were binned 
to 65 A (25 px). The two horizontal segments mark the 
weighted average values within the wavelength range (4500- 
7500 A). The dashed curves trace a Serkowski-Whittet law 
with P mflA =0.15%, 0=169.7 degrees, and 1^=4500 A. The 
upper panel presents the signal-to-noise ratio per resolution el- 
ement in the intensity spectrum. 



Sect. 4. 1 1. This indicates that the Milky Way contribution may 
dominate over that arising within the host. 

4.2. Intrinsic polarization of SN 2005ke 

In order to derive the fundamental parameters of the domi- 
nant axis and their statistical significance, we have performed 
a weighted linear squares fit to the data (accounting for errors 
in both Q and U) in the wavelength range 4000-7500 A, using 
the linear relation U = a + /3 Q. The results are presented in 
Table [3] where we included x 2 P er degree of freedom (DOF), 
the Pearsons correlation coefficient (r rv ), the rms deviation of 
data points from the best fit axis (<x), and the dominant axis 
position angle (8^). 

In addition to the least squares analysis, we have performed 
a principal components analysis. This was achieved construct- 
ing the weighted covariance matrix and extracting its eigenvec- 
tors and eigenvalues, as described in Maund et al. (2010a). The 
parameters retrieved through this analysis are the position an- 
gle of the dominant axis (6></) and the ratio (b/a) of the degrees 
of polarization along the dominant and perpendicular axes. The 
results are presented in Table[3] which includes the centroid co- 
ordinates «e>, (u)). 

Following Wang et al. (2003), we computed the compo- 
nents of the ISP-subtracted polarization along (Pj) and perpen- 
dicular (P ) to the dominant axis, which are obtained by rotat- 



ing counterclockwise the observed Q - U coordinates system 
by the angle that defines the dominant axis (106.1 degrees). 
The result is presented in Fig. [7] where we have also plotted the 
main line identifications, following Taubenberger et al. ([2008). 

The first fact that emerges from this decomposition is that 
the bulk of the polarization is aligned along the dominant axis. 
The orthogonal component P is statistically consistent with 
a null value across the whole wavelength range. One possible 
exception is seen at the wavelength corresponding to the ab- 
sorption trough of the 1 line at about 7570 A. This deviation 
slightly exceeds the 3-sigma level (P a — 0.50±0.15%), and it 
takes place in a region which is partially affected by a strong 
telluric feature (see the telluric absorption spectrum in Fig. [7] 
upper panel), and therefore it needs to be considered with cau- 
tion. Taken at face value, this would imply that oxygen is po- 
larized along a direction orthogonal to that displayed by the 
underlying continuum. 

The second remarkable property revealed in Fig. [7] is the 
relatively low polarizations corresponding to the main absorp- 
tion features such as Si 11 /16355 and the NIR Ca n triplet, with 
respect to that typically observed in Type la SN a week be- 
fore maximum light (Wang & Wheeler 2008). Although there 
are some localized, weak peaks in Pj (some of which are sta- 
tistically significant), the polarization is always below ~0.5%. 
The polarization is dominated by the continuum, which shows 
an increase from blue to red peaking at about 7200 A. This 
is reminiscent of what was found for SN 1999by (Howell et 
al. 120011 1. and establishes a clear link between these two sub- 
luminous events (see Fig. [8] top and middle panels). In addi- 
tion, the continuum polarization marks a clear difference with 
respect to core-normal events. In core-normal events, the con- 
tinuum polarization in the range 6400-7000 A is always below 
0.2% (see for instance SN 2006X, Patat et al. l2009al l. The over- 
all differences between the subluminous events (SN 1999by, 
SN 2005ke) and the core-normal SN 2006X are illustrated in 
Fig. [8] which presents the polarization spectra along the re- 
spective dominant axes. 

The most marked feature in the dominant-axis component 
is the peak corresponding to the Sin /16355 absorption. The 
polarization profile, decomposed along and orthogonally to 
the dominant axis is shown in Fig. [9] The polarization peak 
(P d ~0.39±0.08%) occurs at -13,300 km s -1 , while the mini- 
mum of the absorption trough is at a lower velocity (-10,800 
km s~'), as is typical of Type la SNe (see for instance the well 
studied case of SN 2006X; Patat et al. |2009al >. The polarization 
angle across the line is constant (within the errors), and equal to 
9d, suggesting that the line forming region and the photosphere 
share the same geometrical properties. The signal-to-noise ratio 
is insufficient to tell whether the polarization angle is constant 
across the NIR Ca 11 triplet as well. 

5. Models for SN 2005ke 

In this section, we compare the observations of SN 2005ke 
with possible theoretical scenarios. As in previous studies, 
the results have been obtained using our hydrodynamical ra- 
diation code, HYDRA. This allows us to solve the hydrody- 
namics including nuclear reaction networks for the explosion 
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Fig. 8. Comparison between the dominant axis polarization 
spectra of SNe 1999by (top; Howell et al. |20Qll , 2005ke (mid- 
dle), and 2006X (bottom; Patat et al. |2009aD . The dominant 
component for SN 1999by was computed using #</=80 degrees 
(Howell et al. 2001 ). For presentation, the original 1999by data 
were binned in ~40 A wide bins (11 pixels). The dashed hor- 
izontal line marks the zero-polarization level. The thin (blue) 
lines trace the observed flux spectra. No reddening correction 
was applied to SN 2006X. 



and radiation transport, including detailed atomic models for 
light curves, flux and polarization spectra. The methods are de- 
scribed in Hoflich ( fFWTl [19951 Howell et al. d200TI Hoflich 
et al. (120021 [20061 |2009l , Wang et al. d2006bl) . For calculat- 
ing the polarization spectra in our Monte Carlo module, we use 
the given density, temperature and chemical structures includ- 
ing the non-LTE occupation of the atomic models. We can do 
so because polarization will not affect the rate equations within 
the accuracy of our approximations. In SNe, most of the intrin- 
sic polarization is due to Thomson scattering; nevertheless, our 
calculations take into account also Rayleigh scattering in lines, 
under the assumption of complete redistribution. Thus, the 
scattering matrix is approximated by a linear combination of 
Rayleigh and isotropic shift in phase (Hamilton 1 19471 Domke 
& Hubeny 119881 1 with the size of polarization given by the 
angular momentum of the lower atomic state (Chandrasekhar 
[19601 . 

In the first part of this section, we will discuss the ori- 
gin of the polarization detected in SN 2005ke. This is done 
in the framework of a subluminous delayed-detonation (DD) 
model of a Chandrasekhar mass White Dwarf (WD) because 
many properties of SN1991bg-like events can be understood 
within this family of models. We note that the deflagration- 
detonation transition may occur both during the initial expan- 
sion phase or after a pulsation (Khokhlov, Miiller & Hoflich 
[19931 Hoflich [19951 Hoflich et al. 120021 . We choose a sce- 
nario with MwD-Mch because only these models burn under 
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Fig. 9. Lower panel: polarization profile of Si n /16355 at the 
combined pre-maximum epoch. Solid and dashed lines trace 
Pd and P a , respectively (12 px bins, ~35 A). The dotted-dashed 
curve is the unbinned flux spectrum arbitrarily scaled. Upper 
panel: polarization position angle. The horizontal dashed line 
indicates the dominant axis angle Oj. 



sufficiently high densities to undergo electron capture as in- 
dicated by observations (Hoflich et al. 120041 Motohara et al. 
120061 Fesen et al. [20071 Maeda et al. [20TTb . We note that the 
evidence for electron capture has been found in core-normal 
and moderately subluminous SNe la, SN 1986G-like events, 
but there is a lack of observations of 91bg-like events. Many of 
the arguments presented here are generic and they are also used 
to discuss possible alternative scenarios including mergers and 
edge-lit sub-Chandrasekhar mass models. 

5.1. Delayed Detonation Model 

We start with a spherical model and remap that structure to 
aspherical grids. We can do so because the deviations from 
sphericity are found to be small. The remapping allows us to 
study various effects, such as rotation. These effects may lead 
to oblate, ellipsoidal structures (Eriguchi 1985) and off-center 
Ni distributions, as may be produced by off-center DD (Livne 
[19991 Gamezo, Khokhlov & Oran [20051 Hoflich & Khokhlov 
120061 R6pke[2007l Kasen et al. l2009l . 

Within the DD scenario, the free parameters and prescrip- 
tions are: 1) the chemical structure of the exploding WD, 2) 
the central density, p c , at the time of the explosion, 3) the de- 
scription of the deflagration front, and 4) the density, p lr , at 
which the transition from deflagration to detonation occurs. In 
this work, the spherical structure of the exploding C-O WD is 
based on a solar metallicity model star with 5 M0 at the main 
sequence. Through accretion, this core has been grown close to 
the Chandrasekhar limit (see model 5p0y23z22 in Dominguez 
et al. 2001). At the time of the explosion of the WD, its central 
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Fig. 10. Density (dotted) and velocity (solid) as a function of 
enclosed mass (upper panel), and abundances of stable isotopes 
(lower panel) as a function of the expansion velocity for DD- 
model 5p0z22.8, with a transition density of p, r = 8 xlO 6 g 
cm -3 . In addition, 56 Ni is given (from Hoflich et al. 2002). 



density p c is 2.0xl0 9 g cirT 3 and its mass is close to 1.37M . 
The transition density p tr has been identified as the main factor 
that determines the 56 Ni production, and thus the brightness of 
a SNe la (Hoflich [19951 Hoflich, Khokhlov & Wheeler [J9951 
Iwamoto et al. 1999). The deflagration-to-detonation transition 
density is p, r =8xl0 6 g cm~ 3 . This value, which is lower than 
that typical of core-normal models (2xl0 7 g cm 4 ), is required 
to produce a subluminous explosion (Hoflich et al. 2002). We 
employ this model (5p0z22.8, Hoflich et al. 2002), because it 
reproduces the light curves and spectra of SN 1999by from 
the optical to the NIR. We note that, although some 0.3M Q are 
burned during the slow deflagration phase, only ~ 0.09M o of 
56 Ni are produced because of electron capture. The light curve 
reaches its maximum about 14.6 days after the explosion, at 
M v ~ -17.2. 

The final density, velocity and chemical structures are given 
The density and velocity profiles are very similar 



in Fig. 10 



to core-normal SNe la, because the specific nuclear energy is 
nearly the same when burning from C-O to Si- or Fe-group ele- 
ments; however, more partially burned or unburned material is 
left in the subluminous models. Due to the larger pre-expansion 
during the deflagration phase, most of the star is burned to Si- 
group elements rather than to 56 Ni. This structure is consistent 
with those obtained by inverse spectral analysis (Mazzali et al. 
119971 ). A further feature that distinguishes subluminous mod- 
els from core-normal DD models is the extended region of un- 
burned C-O. As discussed earlier (see also Fig. [9]), the Si line 
profile extends in the blue up to about 16,000 km s _1 , consistent 
with the chemical structure of this model. The photosphere is 



at a velocity of about 10,000 km s 1 , i.e. well within the Si-rich 
zone. 

5.2. Flux spectra 

In Fig. [TT] we compare spectra of SN 2005ke at 8 days be- 
fore maximum with theoretical spectra at about day 9 after the 
explosion assuming an ellipsoidal geometry with an axis ratio 
b/a=0.85. 

The model is taken directly from the time evolution of the 
initial explosion model without further tuning. The adopted b/a 
axis ratio allows us to reproduce the size of the polarization in 
the red for inclinations between 30 and 60 degrees. 

Overall, the theoretical spectra provide a good match to the 
observed flux spectra of SN 2005ke: the overall energy distri- 
bution, the presence and Doppler shifts of the lines are well re- 
produced. The spectrum is dominated by single ionized stages 
of Si, S, Ti and the iron group elements. Given the simplify- 
ing assumptions made, one cannot expect a perfect fit. For one, 
we assumed that the time evolution is identical to the spherical 
case, and we calculated the asymmetric profiles for a specific 
point in time. We had to do so because detailed atomic mod- 
els for the ions are beyond our current computational means. 
This, however, is a good first-order approximation, again be- 
cause the deviations from sphericity are small. The theoretical 
line strengths are slightly larger for most inclinations, which 
may indicate a larger than predicted population of excited lev- 
els of iron-group elements. In addition, the models have more 
line blocking at about 5200 A, and the Tin feature at about 
4800 A is weaker than observed, as is the Can emission at 
about 3800 A. From previous systematic studies of the temper- 
ature effect on spectra (Nugent et al. 1995), we estimate that 
the local temperature in our model is too high by about 1000 
K. 



5.3. Polarization spectra for oblate geometries 

The overall polarization agrees with the data of SN 2005ke 
(Fig. [TT] lower panel). In the optical, P increases with increas- 
ing wavelength, and shows a maximum around 7000 A. The 
largest discrepancy between the models and the data is seen at 
about 5200 A, where the calculations predict little or no polar- 
ization, while the observations show P « 0.2 to 0.3%, i.e. well 
outside the error bars. A missing feature is the depolarization at 
about 8000 A which may be attributed to shell interaction with 
the environment (Gerardy et al. 2004 Quimby et al. 2006b), 
but has not been included in this model. To fit such a feature, 
two additional parameters would have to be added, the orien- 
tation and the mass of the circumstellar matter. Note that the 
seemingly high observed polarization below 3600 A is affected 
by very large uncertainties due to the low flux levels. 

In the following, we will discuss the polarization results. 
In general, polarization provides information about aspheric- 
ity on scales comparable to the mean free path of photons. 
Polarization is decreased if the radiation field becomes more 
isotropic, e.g. by thermalization (Hoflich Tl 99 II ). 
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Fig. 11. Model flux (top) and polarization (bottom) spectra of the subluminous DD model 5p0z22.8 at 7 days before maximum 
light at three inclination angles in comparison with the data of SN 2005ke (~8 days before maximum). We assumed an aspher- 
ical density structure with an axis ratio b/a=0.82. We show the observations (black, with error bars for the polarization). For 
presentation, the observed polarization data were binned to 78 A (30 pix). 



The origin of line and continuum polarization in rapidly ex- 
panding atmospheres dominated by Thomson scattering can be 
represented by three idealized scenarios (H6flich |199U Hoflich 
119951 Kasen et al. I2003I I. These are: 1) aspherical density dis- 
tributions e.g. by rapid rotation of a single WD or binary WD 
mergers, 2) anisotropy in the radiation field, e.g. aspherical 
56 Ni distributions, or 3) selective line blocking of the under- 
lying "photo-disk", e.g. by aspherical chemical distributions as 
caused by off-center DDs. We first consider the latter mecha- 
nism. The alternative polarization mechanism will be consid- 
ered in S ections [6. 1 1 and I6T21 

5.3.1 . Line polarization 

Selective line blocking requires a non-radial chemical gradi- 
ent in the line forming region and can occur even though the 
Thomson scattering photosphere may be spherical. As a re- 
sult, the "photo-disk" shows local polarization increasing with 
the distance from the center, because perpendicularly scattered 



photons and forward scattered photons are 100% and 0% polar- 
ized, respectively (note that, at the photosphere, not all photons 
are traveling radially even at the rim of a disk). For a sphere, the 
components of the polarization field cancel out. If the chemical 
distribution is asymmetric, the disk may be only partially ob- 
scured at the wavelength of the line. If seen off-axis, the line 
may be optically thick in one part of the disk and optically 
thin in another. Thus, the polarization vectors of the underly- 
ing Thomson photosphere incompletely cancel yielding a high 
line polarization. This effect can cause a line polarization of the 
order of 2% if the line-forming region is detached from the pho- 
tosphere, or for strong lines in rapidly expanding photospheres 
ofSNe. 

Strong polarization in lines is produced by '"asymmet- 
ric covering" of the underlying Thomson scattering-dominated 
photosphere by asymmetries in the abundance distribution. For 
strong lines, the amount of material required to selectively 
block the photosphere may be small. Line polarization is a tran- 
sient phenomenon, dictated by the position of the receding pho- 
tosphere with respect to the boundaries of chemical asymme- 
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tries (see also Hoflich et al. 2006 and the discussion in Patat et 
al. [2009a ). High line polarization levels are commonly seen in 
Ca ii, Mg ii, S ii and Si n in many core-normal SNe at early times 
(Wang et al. l2003l Kasen et al. l2003l Wang et al. l2004ll2006bl 
Patat et al. I2009al Maund et al. 1201 Oat . In these objects, the 
effect is large when the photosphere is close to the chemical 
boundary, e.g. between the Mg/Ne/O and Si/S-rich layers. In 
subluminous events, the corresponding lines block the entire 
photosphere, because the Thomson photosphere is smaller, and 
the lines are formed in an extended region of Si/S (see below). 

Although the models included Thomson scattering on free 
electrons and Rayleigh scattering in lines, the former is larger 
by a factor of 10 and, for the sake of the argument, we as- 
sume that absorption in a line depolarizes the radiation field. 
In SN 2005ke and at -8 days, the lines are formed at about 
10,000-12,000 km s _1 , i.e. far from chemical boundaries (see 



Fig. 10 1. Therefore, line polarization by selective line blocking 
is small for strong lines such as Mg a, Si n and S n. Intrinsic 
line polarization due to Rayleigh scattering is small for most 
line transitions; however, high polarization is sometimes ob- 
served in strong lines if they form in regions of asymmetric 
chemical distributions. The resulting asymmetric blocking of 
the light from the scattering photosphere can result in a high 
net polarization. In the model for SN 2005ke, the large exten- 
sion in velocity space of zones of C burning and incomplete 
O burning causes blocking of the entire photosphere, so strong 
lines forming in that region do not display polarization. Weak 
lines form close to the photosphere and, thus, always result in 
small polarization. In conclusion, the polarization spectrum is 
dominated by depolarization in strong lines. 

5.3.2. Continuum polarization 

In addition to the polarization for line features, it is also impor- 
tant to understand the polarization behavior of the continuum. 
In SNe la, a pseudo-continuum is formed by blending of a large 
number of lines which are mostly responsible for thermaliza- 
tion (Kaip [T977i l. The wavelength dependence of the pseudo- 
continuum is one of the keys to understand the wavelength de- 
pendence of P. The ratio R between true absorption k and total 
extinction, k + o~Thomsan> determines the thermalization optical 
depth, t w y/T/~R and, more important, the ratio between po- 
larizing, i.e. Thomson, and non-polarizing extinction at the ra- 
dius of last interaction, T e . vf /„ e „ „ ~ 1. In Fig. 12 R is given for 
Si and Fe-rich mixtures for conditions typical of Type la SN 
photospheres before maximum light. Photons thermalize and 
are depolarized at Thomson optical depths between 0.1 to 3 
and 0.1 to 1, respectively. The ratio depends strongly on wave- 
length. In the UV, where R is close to unity, the opacities are 
line dominated whereas at longer wavelengths,Thomson scat- 
tering becomes more important. 

A second key for the understanding of continuum polar- 
ization in scattering-dominated atmospheres is the behavior of 
polarization as a function of optical depth and degree of asym- 
metry (Hoflich [199 11 1. In Fig. 13 we show the continuum po- 



larization P(t) for various configurations. Let us first consider 
an oblate ellipsoid. If the absorptive optical depth is small, P 
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Fig. 12. Comparison of the ratio R of line to total opacity for 
conditions typical for the photosphere of a SN la at the ther- 
malization optical depth which corresponds to an optical depth 
of «2 to 3 in Thomson scattering for photons in the optical and 
near IR. Line opacity is given at 10,000 K and for a density 
of 10~ 12 g cm" 3 for the Si-rich (top) and Fe-rich (bottom) com- 
position (from Hoflich 1993), corresponding to the thermaliza- 
tion and last scattering layers in subluminous and core-normal 
models, respectively. In the red (at x 7000A) the average R is 
higher by about 50% in the latter. The scattering optical depth 
at the last scattering radius is about 0.5 and 0.3 for sublumi- 
nous and core-normal SNe la, respectively. As a consequence, 
continuum polarization will be smaller for the same degree of 
asphericity (see text). The dashed lines are the flux distribu- 
tions for Planck functions of the corresponding effective tem- 
perature. 



increases linearly with r JC because of the increasing likelihood 
that a photon undergoes scattering. It reaches a maximum at 
t ~ 1.2 and, then, declines by about 60%, as shown by the 



solid curve in Fig. 13 (upper panel). This decrease in polariza- 



tion at larger optical depth is because, with multiple scatter- 
ing, the radiation field becomes more isotropic. Note that for 
an oblate geometry, P varies with optical depth, but does not 
change sign, i.e. the polarization angle is constant. 

This has two important implications for our understanding 
of polarization spectra: 1) the continuum polarization degree is 
not a linear function of the opacity in the pseudo-continuum, 
and 2) the continuum polarization angle may change with 
wavelength even for axially symmetric configurations (i.e. P 
may be not positive definite) for prolate structures or off-center 
energy sources (see Sect. 



6.1 



For oblate ellipsoids the wave- 
length dependencies of opacity and continuum polarization are 
related. Therefore, the wavelength dependence of P seen in 



SN 2005ke (Fig. 1 1 1 can be understood in terms of the vari- 
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Fig. 13. Continuum polarization P for aspherical configura- 
tions in scattering dominated atmospheres with an electron dis- 
tribution oc r~ 3 . Upper panel: maximum linear polarization as a 
function of the optical depth for oblate (solid blue) and prolate 
(dashed red) ellipsoids. Lower panel: P as a function of the in- 
clination for a source located off-center, at a distance D-R p /,/R 
of 3 (solid blue) and 1.5 (dotted red) (adopted from Hoflich 
et al. [19951 Note that D = 3 corresponds to SN 2005ke at 
day -8, when the photosphere R p h is formed at about 10,000 
km s _1 , and assuming the maximum off-center 56 Ni distribu- 
tion suggested by Maeda et al. d20 1 1 1 > for core-normal SNe la. 



ation of the optical depth with wavelength of the photosphere 
for a Si-rich chemistry, which is relevant to subluminous SNe la 
(as opposed to Fe-rich, core-normal events). At short wave- 
lengths, lines form a pseudo-continuum and dominate the opac- 
ity. As a consequence, P is low. Going to longer wavelengths, 
on averag^]/? decreases to attain a minimum at about 7000 A 
(R w 0.5) and P reaches its maximum. 

Between 5000 and 5400 A, blending depends sensitively on 
the temperature, as already discussed in the context of the flux 



spectrum (see Sect. 5.2 1. Therefore, the discrepancy between 
the model and the data seen around 5200 A (see Fig. 1 1 lower 
panel) may be explained in terms of temperature changes, be- 
cause the energy differences between levels in heavy elements 
is comparable to the thermal energy. 

In general, for oblate geometries, the polarization spectrum 
of SN 2005ke can be explained in terms of an intrinsically po- 
larized continuum depolarized by lines. 



3 The curves shown in Fig. 12 1 are a realization at a single value of 
velocity, and do not include the blurring produced by the velocity gra- 
dient within the decoupling region. Although R is sensitive to wave- 
length on small wavelength scales in Fig. [T2| l, the photosphere forms 
over a range of velocities, leading to a significant velocity smearing in 
the effective value of R in a given wavelength range. 



6. Discussion 

SN 2005ke is a close spectroscopic and photometric twin of 
SN 1999by, and both show similar polarization properties, with 
respect to the amplitude of the polarization (~0.6%) as well 
as wavelength dependency. SNe 1999by and 2005ke have very 
distinct polarization properties compared to core-normal events 
(see also Fig. [8}, which place them in a fairly isolated region of 
the Am^ - P plot (Fig. 14 1. They also appear as outliers from 



the relation between the Si n velocity gradient and the associ- 
ated polarization (Maund et al. 2010b). 

The difference in appearance of subluminous events is 
caused by low photospheric temperatures in combination with 
layers of unburned carbon and more massive layers of the prod- 
ucts of explosive carbon and oxygen burning. 

The polarization of SN 2005ke can be understood in the 
framework of an oblate ellipsoidal geometry, with a degree 
of asymmetry of about 15%. If this is typical of subluminous 
Type la events, then a spread in the apparent peak brightness 
by about 0.2 mag is expected, just because of the directional 
dependence of the luminosity. 

One important question to ask is what causes the overall 
asymmetry in these objects, as opposed to what is seen in core- 
normal SNe la explosions like SN 2006X (Patat et al. 2009b 
See also Fig. IS}, and whether this can tell us something about 
the explosion mechanism that leads to these events. 

It is tempting to attribute the asymmetry to the initial 
WD rotation, as first proposed by Howell et al. (2001) for 
SN 1999by. As a first approximation, we assume that the 
original degree of asymmetry is conserved during the explo- 
sion. As for mass redistribution, three dimensional calculations 
for SNe la explosions show little coupling between latitudes. 
When starting from a spherical WD, the resulting envelope re- 
mains spherical even if the burning is highly asymmetric (e.g. 
Gamezo et al. 2005). At about 7 days after the explosion, about 
8 days before maximum light, the photosphere is formed in 
a layer enclosing ~0.7 M WD (Fig. [1 0} . Ba sed on the work 
on rotating WDs (Eriguchi & Muller 7l985l [19931 Miiller & 
Eriguchi 1985 ), a rotation velocity close to the break-up veloc- 
ity produces an asphericity of about 15% at layers placed at 
0.7 in mass coordinates. Fast rotation may be reasonable be- 
cause, within the single degenerate scenario, mass is accreted 
from a companion star via an accretion disk rotating at near- 
Keplerian velocities, although this aspect does not differenti- 
ate subluminous SNe from single degenerate models of core- 
normal events. 

If the scenario in which the overall asymmetry is inherited 
from the WD rotation is correct, then one can reverse the argu- 
ment, and speculate that a low (or null) polarization indicates 
a slow (or no) WD rotation. Now, for core-normal SN la, we 
have upper limits to the continuum polarization only before and 
around maximum light. This limit is 0.1 to 0.2% (Wang et al. 
I2006bl Wang & Wheeler [20081 Patat et al. I2009al l. This in- 
dicates that fast WD rotation and SN sub-luminosity may be 
related, marking a clear distinction with respect to core-normal 
explosions. 

At this point one can ask: do we expect a similar degree 
of polarization in normal and subluminous events for the same 
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degree of asymmetry? In core-normal Type la and around max- 
imum light, the thermalization of pseudo-continuum photons in 
the red takes place mostly in layers of the iron group, at vari- 
ance with objects like SN 2005ke, in which thermalisation oc- 
curs in layers of incomplete burning. The ratio R between true 
absorption and total extinction (see Sect. |5.3| l, averaged over 
the region 6750-7250 A, is 0.7 and 0.5 in core-normal and sub- 
luminous supemovae, respectively. There are two reasons for 
the difference in R. Firstly, the thermalization layers changes 



from iron-rich to Si/S-rich compositions (Sect. 5.3 1. Secondly. 



for core-normal SNe la, the thermalization happens at layers of 
larger expansion velocity and, thus, smaller density leading to 
a further increase of R (see Hoeflich et al. 1 19931 their Fig. 1). 

Moreover, larger thermalization in core-normal SNe la re- 
sults in an increased isotropy of the radiation field (Hoflich 
1991). For the same amount of asymmetry, the intrinsic po- 
larization at about 7000 A may be smaller by up to a factor « 2 
for core-normal SNe la (Figs. 12 1, 13 I, but not a factor 5 to 



6 as indicated by early-time spectropolarimetric data. Though 
not conclusive, because of the small number of early-time ob- 
servations for SNe la, there indeed may be a connection be- 
tween rotation and sub-luminosity. In this respect, it is interest- 
ing to note that a relatively massive (1.28+0.05 M Q ) WD, with 
a rapid rotation (P= 13.2s) has been found in a binary system 
(Mereghetti et al. 1201 11 1. The companion star is a mass-loosing, 
hot sub-dwarf, from which the WD is currently accreting ma- 
terial. 

6.1. Prolate geometries and off-center energy sources 

Given the lack of spectropolarimetry time coverage, alterna- 
tive geometries within Met, explosions need to be consid- 
ered, namely prolate ellipsoids and off-center energy sources. 
These can be understood within the same context discussed in 



Sect. 5.3 In prolate geometries, tangential rays dominate the 



polarization for small optical depths causing negative P. With 
increasing optical depth, multiple scattering causes the radia- 
tion field to become more isotropic. This reduces the tangen- 
tial flux and causes P to change sign at greater optical depth 
(Fig{!3] and Hoflich [T99T1H995l . 

In the light of these considerations, we first note that overall 
prolate structures can be ruled out for SN 2005ke. The observed 
polarization is positive at all wavelengths whereas, in the case 
of a prolate geometry, it should change its sign (i.e. angle) at 
wavelengths dominated by Thomson scattering (Fig. |T3]>. 

As for off-center energy sources (i.e. 56 Ni) in Type la SNe, 
these have been inferred from the NIR line profile at 1 .65 fim 
(Maeda et al. 1201 lb . and from observations of the remnant 
of S-Andromedae, which has been attributed to the class of 
SN 1986G-like objects (Fesen et al 120071 . The amount of "off- 
centerness" is less than about 2000 km s~ 1 . The observations of 
remnants and NIR line asymmetries are limited to a few cases 
only, none of which is a SN 1991bg-like event. For DD models, 
a week before maximum, the optical depth of the central region 
is about 50 (Hoflich Tl 9951 ). At -8 days, the inner layers cannot 
contribute to the emission because their diffusion time scales 
are larger than the expansion time. Therefore, it is not possible 
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Fig. 14. Degree of polarization across the Si n /16355 line as a 
function of light curve decline rate. Data are from Wang et al. 
(120071 ) and Patat et al. d2009al ). SN 2005ke is marked by the 
large dot close to SN 1999by. The solid line is a best fit to all 
data, with the exception of SN 1999by, 2004dt, and 2005ke. 
The dotted lines trace the 1-cr level of the intrinsic polariza- 
tion distribution generated by the Monte-Carlo simulation dis- 
cussed by Wang et al. (2007). The shaded area indicates the 
rms deviation of the data points from the best fit relation. 



to distinguish global density asymmetries and 56 Ni asymme- 
tries based on pre-maximum spectra. 

Nevertheless, we regard off-center explosion an unlikely 
source of the polarization observed in SN 2005ke. If an off- 
center explosion is causing the early continuum polarization, 
the polarization should increase until the photosphere has re- 
ceded to the Ni layers at about 2 weeks after the explosion (i.e. 
around maximum light) and, then, the polarization angle should 
flip. Our late time spectropolarimetry of SN 2005ke does not 
show this effect, lending support to the conclusion that global 
density and not energy input asymmetries are at the origin of 
the continuum polarization displayed by this object. 

Chemical asymmetries would cause large polarization in 
lines and no continuum polarization. Therefore, in the case of 
subluminous events, where the lines form far from chemical 
boundaries, the observed continuum polarization requires the 



density distribution to be asymmetric (see Sects. 5.3.1 5.3.2 1 



6.2. Sub-Chandrasekhar Mass Explosions and 
Mergers 

It is still under debate whether core-normal and subluminous 
Type la are a single class (Hoflich et al. 2002 ), or form separate 
groups (Hillebrandt et al. l20T0l Sim et al. I20TT1 Ropke et al. 
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1201 p . Both sub-Chandrasekhar mass WDs and WD mergers 
have been suggested as possible mechanisms. 

In sub-Chandrasekhar mass models the explosion is trig- 
gered by a detonation in the He-layers on top of a low mass 
WD which triggers the detonation of the C-O core. Explosive 
He burning is much more energetic, and the resulting structures 
of all models show either a Ni layer (Hoflich & Khokhlov 1996 1 
Nomoto et al. [19971 Hoflich [19971 Woosley [19971 Woosley & 
Kasen l20TT1 l or, with some fine-tuning, strong Fe- and Ti-rich 
layers (Sim et al. 1201 II Woosley et al. 201 1). As a generic fea- 
ture, this class of models show a 56 Ni core, surrounded by a 
layer enriched in Si-S, subsequent layers of O-Mg-Ne, C-O, 
Ni-Fe-Ti, and an outer layer of high velocity He. Woosley & 
Kasen (1201 11 1 show that maximum light spectra look reason- 
able because the spectra are formed in Si/S layers at a depth 
of about 0.1 M Q . Successful models must, however, agree at 
all phases. In SN 2005ke at day -8, the spectra are formed in 
the very outer layers, but we do not detect any of the chemi- 
cal signatures predicted by any of the sub-Chandrasekhar mass 
models. 

Merger of two WDs (Iben & Tutukov [l984l Webbink [l984l 
Benz et al. [19901 Pakmor et al. [20TU1 [20TTb has recently drawn 
attention because of perceived advantages for solving the Type 
la progenitor problem (Ruiter et al. 12009 1 Schaefer & Pagnotta 
al. I20121 >. There is evidence that mergers may form a subclass 
among SNe la (Hoflich & Khokhlov [T9961 Quimby, Hoflich & 
Wheeler|2002l Pakmor et al. l20T0ll20Tll . 

Merger models can show a layered chemical structure simi- 
lar to SD models of M C / : explosions (details depend on whether 
or not material from the accreted WD is assumed to detonate). 
There are two aspects that qualitatively distinguish merger 
models from SD models: 1. late-time IR spectra and the rem- 
nant to S-Andromeda show flat topped profiles that seem to 
require central densities beyond 10 9 g cirT 3 in order to have 
electron capture time-scales comparable to hydrodynamical 
timescales. This condition is required to have a 56 Ni deple- 
tion in the center (Hoflich et al. 120041 Maeda et al. 1201 11 1: 2. 
merger models tend to produce a significant amount of un- 
burned carbon in the outer layers. Core-normal SNIa show lit- 
tle or no unburned carbon (Marion et al. 2006; Parrent et al. 
1201 II Silverman & Filippenko 120121) . These arguments argue 
against a merger origin for the majority of core-normal (and 
SN1986g-like in the case of S And) SNe la. For the sublumi- 
nous SN 2005ke, these discriminators against merger models 
cannot be applied. Even explosions of single-degenerate WDs 
such as DD-models for SN 1999by show strong carbon lines in 
the NIR down to about 14,000 km s . The existence of carbon 
does not sway the argument one way or the other. Moreover, we 
have no late-time NIR spectra for SN 2005ke or for any other 
subluminous SNe la and therefore no constraint on the central 
density in the progenitor WDs. Although we need a WD with a 
mass close to M c/l for SNe 1991bg, 1999by and 2005ke to ex- 
plain their light curves and spectral evolution (Hoflich 2002), 
we cannot rule out an object with similar but slightly lower 
mass (Mazzali et al. 119971 Taubenberger et al. 12008b . In the 
context of the current work, the relatively high degree of con- 
tinuum polarization observed in the subluminous SNe 1999by 
and 2005ke could imply a high degree of rotation in the pro- 



genitor and hence may be consistent with a double-degenerate 
merger, for which the explosion would be intrinsically asym- 
metric (e.g. Pakmor et al. 1201 11 1. 

We may expect that mergers induce strong differential ro- 
tation and, thus, larger asymmetries well after maximum light 
(Eriguchi & Muller [T985b . However, hydrodynamical models 
for mergers evolved to the phase of free expansion are still 
missing. In general, we would expect large asphericity in the 
central region for mergers, but almost spherical cores for ro- 
tating WDs. In light of these considerations and the available 
data, we cannot rule out a merger for SN 2005ke. 

Time sequences of spectropolarimetric data and late-time 
NIR spectroscopy of subluminous events may provide a more 
definitive conclusion. 

7. Conclusions 

We presented and discussed the pre-maximum spectropo- 
larimetry of the subluminous Type la event SN 2005ke. The 
main observational results can be summarized as follows: 

1 . SN 2005ke shows similar properties to the only other sub- 
luminous event for which spectropolarimetry is available 
(SN 1999by). They differ substantially from core-normal 
SNe la. 

2. Along the dominant axis the two subluminous events dis- 
play a significant continuum polarization (0.6-0.7% at 
peak) that increases steadily toward the red, at variance 
with core-normal events, which show a low continuum po- 
larization (<0.2 %) that is essentially constant in the optical 
band. 

3. The relatively strong polarization associated with absorp- 
tion lines in core-normal events is absent in the two sublu- 
minous objects. 

We interpreted these findings using our hydrodynamical ra- 
diation transfer modeling. The following conclusions emerged 
from this analysis: 

1. The flux and polarization spectra of SN 2005ke are well 
reproduced by an oblate ellipsoidal geometry within a sub- 
luminous delay ed-detonation scenario. 

2. The differences with respect to core-normal SNe la are 
caused by low photospheric temperatures in combination 
with layers of unbumed C and more massive of the prod- 
ucts of explosive C and O burning. 

3. The comparatively large continuum polarization is ex- 
plained in terms of a global asymmetry (~15%), which is 
not present in core-normal explosions 

4. In the two subluminous events, the lines of intermediate 
mass elements form far from chemical boundaries and over 
a large velocity range compared to core-normal SNe la. 
This causes a blocking of the entire photosphere, resulting 
in weak line polarization. 

5. The overall asphericity characterizing subluminous Type 
la may be produced either by a fast WD rotation, or by a 
double-degenerate merger. 
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